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Abstract 

The conductance of a waveguide containing finite number of periodically 
placed identical point-like impurities is investigated. It has been calculated 
as a function of both the impurity strength and the number of impurities using 
the Landauer-Biittiker formula. In the case of few impurities the conductance 
is proportional to the number of the open channels N of the empty waveguide 
and shows a regular staircase like behavior with step heights ~ 2e 2 /h. For 
large number of impurities the influence of the band structure of the infinite 
periodic chain can be observed and the conductance is approximately the 
number of energy bands (smaller than N) times the universal constant 2e 2 /h. 
This lower value is reached exponentially with increasing number of impu- 
rities. As the strength of the impurity is increased the system passes from 
integrable to quantum-chaotic. The conductance, in units of 2e 2 //i, changes 
from N corresponding to the empty waveguide to ~ N/2 corresponding to 
chaotic or disordered system. It turnes out, that the conductance can be ex- 
pressed as (1 — c/2)N where the parameter < c < 1 measures the chaoticity 



of the system. 

In recent years transport in a wide variety of mesoscopic systems has been investigated. 
One of the interesting questions is the behavior of the conductance of disordered identical 
blocks coupled in chain. This problem may arise whenever identical chaotic or disordered 
nanostructures are organized in a sequentially built structure. So far mostly chains with 
completely disordered/chaotic blocks have been studied Our aim in this letter is to 

study a system in which chaoticity of the blocks can be tuned by varying a single parameter 
c ranging from zero (completely regular) to one (completely chaotic). Well known examples 
of such systems are provided by point-like impurities with variable strength placed in regular 
waveguides H||. We study here the low temperature conductance of a waveguide containing 
finite number of periodically placed identical point-like impurities. In this model the con- 
ductance depends on the strength and on the number of impurities. It will be shown that 
when passing from the regular toward the quantum-chaotic case by increasing the impurity 
strength the average conductance (in units of 2e 2 /h) behaves like ~ aN, where a = 1 — c/2, 
interpolating between the conductance of the regular (empty) waveguide N and the uni- 
versal conductance of chaotic systems ~ N/2 determined by Random Matrix Theory |3j 
(RMT). The direct universal link between c, characterizing the chaoticity, and the average 
conductance is our main result. 

We consider an ideal 2D waveguide of width W which is divided into L blocks of length 
a. We place an impurity with Dirac-delta potential U(r) = X5(r — r ) in each block, where 
r is its position within a block (see Fig. 1) and A is the strength of the potential. Inside 
the block the potential is assumed to be zero and the wave functions should fulfill Dirichlet 
boundary condition (■?/> = 0) on the walls of the waveguide. The system is adiabatically 
matched to 2D electron reservoirs on both ends. 

The conductance is given by the Landauer-Biittiker formula 

G{E F ) = ^T, (1) 

where T is the total transmission coefficient of the system at Fermi energy Ep. The total 

2 



transmission can be expressed in terms of the partial transmission amplitudes t nm of open 
modes between the entrance and the exit of the waveguide 

T = ^ S \tnm\ i (2) 
n,m 

where the transmission amplitudes can be calculated from the retarded Green-function of 
the system 

t nm = 2i{k n k m ) l l 2 f dydy'$ n (y)G + (x,y\x', y')§ m {y')- (3) 



Here & n (y) = J sm(irny /W) is the transverse part of the wavefunction of the empty 
waveguide 



^ n {x,y) = <$> n {y)e iknX , (4) 



where k n = — (him/W) 2 for mode index n, y is the transverse coordinate and 

longitudinal coordinates x and x' lie anywhere on the entrance (left hand side in Fig. 1) and 
exit sides, respectively. The mode n is open if Ep > (him/W) 2 /2m. 

The retarded Green-function can be calculated recursively |J. Adding a new impurity 
to a system with L impurities changes the Green-function the following way 

n + f | / a n+( | / a , \ G${x, y\x , y )Gt{x , y \x', y') 

Gl +l {x, y\x ,y) = Gl{x,y\x ,y) + A ~— — ■ , (5) 

1 - XGl{x ,y \x ,y ) 

where the strength A should be renormalized 

1 + XGt(x , y \xo, yo, -E F ) ' 

in order to make the calculations with the delta potential well defined 0. The recursion 
starts with the retarded Green-function of the empty waveguide 

G+{x,yW,y\E F ) = ^ y ^ y 'K ^\ (7) 

n=l *1Kn 

where n max — > oo. In numerical calculation n max should be much larger than the number of 
open modes iV = IntikpW jixTi). 



First we study the behavior of the conductance of the system with increasing number of 
blocks L each containing one impurity of fixed strength A. For the empty ideal waveguide 
(L = 0) the transmission as a function of the Fermi wavelength kp = y/2rnEp/h is an ideal 
staircase with steps T = N, where N is the number of open modes. Adding one impurity 
(L = 1) slightly modifies the staircase |1(J and the transmission remains still close to N. By 
adding more impurities the staircase disappears and an irregular structure emerges as it is 
shown in the upper part of Fig. 2 (solid line) for L = 10. The transmission for L — > oo at a 
fixed Fermi wavelength decreases toward an asymptotic value which is the transmission of 
the infinitely long periodic system. 

The eigenmodes of the infinitely long periodic system are the Bloch functions 

tpq,n( X ^y) = e iqX U q , n {x,y), (8) 

where u g>n (x + a,y) = u g>n (x,y) and —n/a < q < 7r/a with eigenenergies e n (q) forming 
a band structure. In this case the transmission is given by the number of Bloch modes 
Nb propagating in positive x direction at Fermi energy Ep, which can be determined by 
counting the number of solutions (n*, <&), i = 1, Nb for qi > of 

e ni {qi) = E F . (9) 

In the lower part of Fig. 2 the band structure for a given value of A is shown for q > 0. Here 
one can count the number of intersections of bands with kp yielding iV^ and the transmission 
T = Nb of the infinite system is plotted on the upper part of Fig. 2 (dashed line). We can 
see that the transmission, even for L = 10, approaches integer values corresponding to the 
transmission of the infinite system. The transmission drops to zero when kp is in a gap of 
the band structure and it jumps whenever kp is on the edge of a band. 

Next we study the L dependence of the average transmission < T >n corresponding to 
a fixed number of open modes N, i.e. for Fermi wavelengths from kp = N to kp = N + 1 
(hereafter measured in units of nh/W), 

rN+l 

<T>at= / T(k F )dk F . (10) 
Jn 
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In Fig. 3 < T >io is shown as a function of the number of blocks L for different values of 
A. The average transmission as the function of L decreases exponentially and can be well 
approximated with the expression 

< T > N = (N - (A, N))e~ L ^ x ^ + T^X, N), (11) 

where £ characterizes the rate of decay and is the number of modes remaining open for 
L — > oo. We can interpret £ as a partial localization length in the sense that the conductance 
decays exponentially with the size of the system, however it does not decay to zero, like in 
case of strong-localization, due to the periodicity of the system. 

In Fig. 4 we plotted the asymptotic transmission as a function of the number of 
open modes N for different values of A. The transmission seems to scale approximately 
linearly with N. For large A it is in good agreement with the prediction of the Circular 
Orthogonal Ensemble (COE) of RMT for large N 

N 1 

Tcoe « T - 4, (12) 



worked out for chaotic cavities with time reversal symmetry in Ref. [II |. For other values 
of A one can fit 

Too(A, N) = a(X)N + f3(\). (13) 

In Fig. 5 the fitted parameters a (A) and /3(A) are shown. One can see that j3 is ap- 
proximately independent of A and takes the value —1/4, the weak localization correction, 
predicted by RMT. This means, that this weak localization correction sets in if back scat- 
tering is possible (A ^ 0), however it does not depend on the details of its mechanism. 

The parameter a decreases from a ~ 1 to a ~ 1/2 with increasing A. The A dependence 
of a can be explained by making a simple assumption about the relation between the con- 
ductance and the chaoticity of the system as we are going to show next. The Hamiltonian 
in systems which are neither completely regular nor chaotic can be divided into a regular 
Hr and a chaotic He part 



H = H R + H C . (14) 

In our case Hr is the Hamiltonian of a block of the empty waveguide and He is associated 
with the delta impurity in it. The regular part can be modeled with a diagonal matrix with 
random elements and mean level spacing A. The chaotic part can be modeled with a random 
Gaussian matrix drawn from the suitable orthogonal ensemble (GOE, GUE etc.) reflecting 
the symmetry of the system. The variance of the elements of the random matrix should 
coincide with the variance of the matrix elements of the chaotic part of the Hamiltonian 



i.e. a 2 = | < i\Hc\j > | 2 . The energy level statistics of such systems is transitional between 
Poissonian and tpure RMT. The chaoticity of the system, reflected in the level statistics, is 
determined by the transition parameter 7] = na/A whose universality was shown in many 
applications [fL^] . The appropriate measure of the chaoticity of a system is the rescaled 
parameter defined in Ref. 



(15) 



ranging from to 1. The value c = corresponds to a completely regular system (the empty 
waveguide) and for c = 1 the system is completely chaotic (A — > oo) in a quantum- chaotic 
sense. Analogously, in our numerical calculation (W = a = 1), the chaoticity parameter c 
can be defined as 

c(A) = A/4 . (16) 
y/1 + (V4) 2 

Using semiclassical considerations for large N we can make connection between the chaoticity 
parameter c and the parameter a characterizing for N —>■ oo. Semiclassically a part of the 
trajectories is bouncing between the walls only, while the rest is scattered by the impurities 
irregularly. Accordingly, the transmission can be approximately partitioned into a regular 
and a chaotic part. We can assume that the effective number of chaotic modes is iV • c while 
the number of regular modes is N ■ (1 — c). For regular modes the propagation is ideal, 
free from reflection and their contribution to the transmission is N ■ (1 — c). On the other 
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hand, according to the RMT result flT2|), only half of the chaotic modes contribute to the 
transmission and we get 

Too & N • (1 — c) + • c, (17) 

yielding 

« = l-c/2. (18) 
In Fig. 5 we compared the numerically fitted values of a(A) (see Eq. (|13|)) to 1 — c(A)/2 



see Eq. ([lq) ) and we have found excellent agreement in the whole range A = 1 — 500 



supporting our new formula (jT8|). 

We conjecture that our new formula ( |I8~D is valid for the transmission of a much broader 
class of systems. It is expected that it works for the transmission through a system, attached 



to leads, whose level statistics (without leads) is transitional |H| in between Poissonian and 



RMT and can be characterized with < c < 1, since the assumption made above are valid 
in this case too. This makes it possible to determine the chaoticity of a system from the 
relation 

c = 2-(l-a), (19) 

by estimating a from the conductance. We hope that further numerical and possible exper- 
imental studies will support the results presented here. 
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FIGURES 



x 

FIG. 1. Schematic picture of the waveguide containing periodically placed impurities. The 
distance between the impurities is a and the width of the guide is W . The entrance is on the left 
hand side. In the numerical calculations we choose W = a = 1 and the position of the impurities 
was xq = 0.65 and yo = 0.84 with respect to the left lower corner of the unit cell. 
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FIG. 2. In the upper part the transmission vs. the Fermi wavelength in units of hir/W is 
plotted for a system with L = 10 blocks (solid line) and for an infinitely long system L = oo 
(dashed line). A = 5 in both cases. The rapid fluctuations of the conductance are on the scale 
Ak ~ 1/La caused by the multiple scattering between the two ends of the device. In the lower 
part the band structure of the infinitely long system is plotted for the same parameter A as above. 
The vertical axis is the Bloch number q and k n (q) = (2me n (q)) 1 ^ 2 /% has been plotted in units of 
hir/W. This representation makes it possible to get the transmission of the infinitely long system 
in the upper part from the bands of the lower part graphically 
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FIG. 3. The average transmission for N = 10 open modes as a function of the number of blocks 
L at different values of impurity strength A and for W = a = 1. < T >io decreases from 10 (the 
transmission of the empty waveguide) exponentially as a function of L for all A. Around L « 50 
it reaches an asymptotic value depending on A. For large A (chaotic case) this asymptotic value is 
close to N/2 = 5 predicted by Random Matrix Theory. 
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FIG. 4. The asymptotic transmission as a function of the number of open modes N of the 
empty waveguide for different values of A in the range of 1 — 500. The solid line represents the 
COE result T = N 2 /(2N + 1). 
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FIG. 5. The parameters a and (3 (see Eq. (|l~3|)) of the linear fit to the data of Fig. 4 for 
different values of A in the range of 1 — 500. a, parameter a vs. A (diamonds) and our prediction 
a = 1 — c(A)/2 (see. Eq. (|i~7|)). The errorbars represent the error of the linear fit. b, parameter (3 
versus A (diamonds) and the weak localization correction —1/4 (solid line). 
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